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SUMMARY 

This report describes a study of the flow field, spray penetration, vapor- 
ization, and combustion In two-stroke diesel engines. Fuel injection begins 
at 345° after top dead center (ATDC), and n-dodecane is used as the liquid 
fuel. Arrhenius kinetics is used to calculate the reaction rate term in the 
quast-global combustion model. When the temperature, fuel, and oxygen mass 
fraction are within suitable flammability limits, combustion begins spontane- 
ously. No spark is necessary to ignite a localized high-temperature region. 
Compression is sufficient to increase the gaseous phase temperature to a point 
where spontaneous chemical reactions occur. Results are described for a swirl 
angle of 22.5°. 


INTRODUCTION 

This report describes the flow field, spray penetration, vaporization and 
mixing, and combustion in a two-stroke diesel engine. This study of the firing 
diesel engine uses the same geometry, initial conditions, and boundary condi- 
tions as were used for the motored case studied previously (ref. 1). Liquid 
spray combustion begins at 345° after top dead center (ATDC) (end of compres- 
sion stroke) when fuel injection begins. The Eulerian conditions of the gas 
phase are identical to those used In the motored case (ref. 1) and are not 
repeated here. The initial conditions of the liquid phase are also the same 
as those used in the motored case but they are repeated in the next section. 
Arrhenius kinetics is used to calculate the reaction rate term. The results 
show the firing cycle to be composed of three periods: (1) ignition delay, 

(2) explosion, and (3) controlled diffusive burning. The ignition delay 
period is characterized by fuel vaporization and extensive mixing in an enve- 
lope around the spray. Explosive combustion begins abruptly throughout the 
partially premixed envelope once the temperature of the engine cylinder rises 
to within the fuel flammability limits. Large pressure gradients and pressure- 
driven flow dominate the combustion and droplet dynamics during the explosion, 
resulting in rapid droplet vaporization and dispersion. Large regions of fuel 
vapor develop near the injector. After a short period of explosive premixed 
burning, more controlled diffusion fuel flames develop. Vortices generated by 
the injection of liquid fuel and by several localized combustion regions par- 
tially mix the segregated regions of fuel and oxygen. During expansion, large- 
scale vortices gradually decay as piston-driven flow develops. Although 
localized combustion continues throughout the expansion stroke, insufficient 
large-scale motion results in incomplete mixing of the fuel and oxygen. 


PROBLEM FORMULATION 


The piston-cylinder configuration analyzed herein consists of a (1) cir- 
cumferential Intake port located at the cylinder wall that Is covered and 
uncovered by the piston; (2) an Infinitesimally thin exhaust valve, which is 
located at the center of the cylinder head and penetrates the cylinder at 
specified crank angles; and (3) an injector located at the cylinder centerline 
and cylinder head (refs. 1 to 3). The geometric characteristics of this 
two-stroke diesel engine are presented in figure 1. 

Descriptions of the mean flow field, the two-equation model of turbu- 
lence, and the discrete droplet model can be found in references 1 to 3. A 
one-step irreversible chemical reaction of the Arrhenius type is used to model 
the conversion of fuel <C] 2 H 26 ) and oxidizer (C> 2 > into combustion products 
(CO 2 and H 2 O). The Arrhenius kinetics oxidation reaction rate was assumed to 
be proportional to a temperature-activated exponential term; the values of the 
pre-exponential factor, the activation energy, the powers of fuel, and oxidizer 
concentrations are taken from reference 4. 


PRESENTATION AND DISCUSSION OF RESULTS 

Table I summarizes the engine variables investigated in this report. The 
firing study uses the same geometry, initial conditions, and boundary condi- 
tions as were used by the motored case. The compression period is identical 
to that of the motored case (refs. 1 to 3) and is not repeated here. The 
study begins at 345° ATDC with the beginning of fuel injection; the final out- 
put of the compression stage is read and used as the initial conditions. 

Three distinct combustion periods exist in the cylinder during one cycle. 
The first period, commonly called the ignition delay period, lasts from 345° 
ATDC until approximately 353° ATDC for the geometry being modeled at 1000 rpm; 
the length, however, changes significantly depending on rpm, geometry, intake 
boost, compression ratio, etc. During this time the liquid fuel vaporizes as 
it did in the motored study; however, combustion is Insignificant because of 
the lack of energy to initiate the chemical reaction. Two effects cause the 
delay - (1) the vaporization cooling of the regions of partially mixed fuel 
and oxygen and (2) the insufficient compression of the cylinder gases. During 
this period significant regions of nearly stoichiometric mixture develop in an 
envelope around the spray region. Gradually compression and heat transfer 
from the warmer regions raise the temperature of these regions giving rise to 
the second period. 

An explosive situation exists for about 2° when the temperature of the 
partially premixed region initiates a rapid combustion. The combustion quickly 
consumes the partially premixed fuel in the envelope around the liquid fuel. 
Rapid pressure and temperature rises occur throughout the cylinder. After all 
the initially short premixed regions have been burned, the combustion relaxes 
into the third period. 

This third combustion period is characterized by diffusion-controlled 
flames. The gaseous fuel is located primarily in the injector region, near 
the injection vortex, where the droplets vaporized. The oxygen is located 
around the fuel-rich region, throughout the rest of the cylinder. Combustion 
proceeds but at a relatively slow rate as the fuel and oxygen diffuse 
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together. A diffusion flame envelopes the fuel-rich center and slow and more 
controlled combustion continues until the blow-down period begins. 

Figures 2 to 4 show the ignition delay period. Figures 2 and 3 show the 
velocity, fuel, oxygen, and carbon dioxide plots at 350° ATDC. Figure 2 shows 
the development of the injection-driven vortex and the scraping vortex along 
the cylinder wall. The major feature in the flow is the momentum exchange 
between the liquid and gaseous phases. Figure 3 shows the localized region of 
gaseous fuel mass. Fuel displaces the oxygen near the injector. Note that no 
appreciable combustion occurred anywhere in the cylinder (fig. 3). 

Figures 4 and 5 show the hydrodynamic features at the last instant before 
the explosion period begins. The velocity in the injection vortex (fig. 4) is 
well developed as the liquid-gas momentum transfer drives the flow. In fig- 
ure 5, the effects of vaporization cooling are seen near the injector. The 
temperature is uniform throughout the cylinder with sharp temperature gradi- 
ents near the walls. 

Figures 6 and 7 show the crankshaft degree during which explosive combus- 
tion is present in the chamber. These figures also show the velocity, tempera- 
ture, pressure, carbon dioxide, oxygen, and fuel at 355° ATDC. Figure 7 shows 
a nearly unidirectional velocity field during the explosion. No injection or 
scraping vortex can be seen as the explosive pressure rise drives all aspects 
of the flow. Figure 7 shows the temperature rises in the cylinder. A large, 
partially premixed region located near the cylinder head at the interface 
between the gaseous fuel and oxygen explodes, causing a rapid temperature rise 
in that region. Note the relatively large region over which temperature has 
increased, indicating a large amount of burning and, also, that a large amount 
of partially premixed fuel and oxygen was present at 354° ATDC. Large pres- 
sure gradients exist between the cylinder wall and the centerline, and the 
global pressure in the cylinder has increased by a factor of three. The car- 
bon dioxide is a major combustion product in the burned region near the injec- 
tor. The fuel air mixture originally associated with the injection vortex was 
premixed to within the flammability limits of the fuel and burned during the 
explosion period. Levels of oxygen near the injector are substantially 
reduced, and most of the gaseous fuel outside the dense liquid spray core 
region is consumed in the explosion. This explains the large region over 
which a substantial temperature rise was observed (fig. 7). The regions near 
the cylinder wall where no combustion is present Increase in temperature due 
to the adiabatic compression associated with the large cylinder pressure rise. 

Figures 8 and 9 show the velocity and temperature profiles, respectively, 
at 356° ATDC. Small effects of the explosion period can still be seen here. 
Figure 8 shows a nearly unidirectional rebound from the cylinder wall toward 
the centerline. A small injection vortex is still observed near the Injector 
region, and recirculation regions are present when the flow strikes the center- 
line as well as in the regions near the cylinder wall. Figure 9 shows that 
further burning near the centerline occurs. 

After a short period of 2° of explosive premixed combustion, the burning 
slows into diffusion-controlled combustion. Figures 10 and 11 show the 
velocity and temperature, respectively, at 357°. Figure 10 shows the develop- 
ment of two large combustion-driven vortices and the redevelopment of the 
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Injection vortex. One vortex is located on either side of an interface line 
dividing the fuel-rich center from the oxygen-rich region of the cylinder. 

The fuel-side vortex is driven by the combustion near the piston face Fig- 
ure 11 shows the region near the piston face that is burning. The remains of 
the large premixed regions near the center of the cylinder are burning and 
drive the two large recirculation zones. 

ocoi 9Ur j S -,rno t0 15 show the veloc1t ?> temperature, and carbon dioxide plots 
at. 358 and 359 . The velocity plots indicate the further development of the 
injection-driven vortex. The temperature profiles indicate burning in the in- 
jection vortex near the centerline and at the interface between the oxygen and 
fuel vortices near the cylinder head. 

Figures 16 to 18 show the velocity, temperature, carbon dioxide, oxygen 
and fuel profiles at 360° ATDC. Figure 16 shows the continuing development of 
the injection-driven vortex. The temperature profile (fig. 17) shows only 
weak combustion near the injection-driven vortex. Figure 18 shows the reason 
for the lack of combustion. The fuel located predominantly near the injection 
vortex Is gradually drawn into the vortex located in the fuel-rich region. 

The oxygen is located primarily near the vortex located in the oxygen-rich 
region. In figure 18, the result of the combustion can be seen near the cylin- 
der head next to the oxygen-rich region where the oxygen is depleted and the 
carbon dioxide is abundant. This explains the high-temperature region that is 
present In the cylinder head near the oxygen-rich region. 

Figures 19 to 22 show the velocity profiles at 1° to 4° ATDC, respec- 
tively. The velocity near the injector develops into a strong injection-driven 
vortex. Both combustion-driven vortices decay substantially, and combustion 
has effectively stopped. 

Figures 23 and 24 show the velocity and temperature at 5° ATDC. The 
velocity profiles show that on the last degree of Injection the liquid 
momentum-driven Injection vortex is again established and is driving the flow 
n the region near the injector. Combustion is nearly extinguished and no 
longer drives the mixing of fuel and oxygen. The temperature profile indi- 
cates oxygen-rich outer regions and fuel-rich center regions. Correspond- 
ingly, the temperature is hottest at the interface and cooler on each side. 

The oxygen-rich region is cooler since no combustion has occurred there. 

There are no large-scale motions that are mixing the fuel with the oxygen. 

The fuel is swept into the center of the cylinder by the injection vortex but 
most never reaches the oxygen-rich region. 

Injection ends at 5° ATDC, and the injection-driven vortex moves more 
toward mid-cylinder. Figures 25 and 26 show the velocity and temperature 
profiles at 10° ATDC. The injection vortex has adjusted itself to the lack of 
droplets near the injector, and there are indication of the piston beginning 
to drive the flow. Both the fuel-rich region and oxygen-rich side vortices 
are stretched by the presence of the injectors vortices. Rapid expansion away 
from the burning causes a separation in the flow (fig. 25). A rise in tempera- 
ture develops near the cylinder head causing new combustion near the cylinder 
head (fig. 26). 

Figures 27 and 28 show the previously mentioned variables at 15° ATDC. 
Figure 27 shows the piston beginning to drive the flow. The combustion 
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present at 10° ATDC has decayed. Figure 28 shows a temperature peak near the 
center of the cylinder where combustion raises the gaseous temperature. Bound- 
ary layers are developing along the three walls. No large-scale motion exists 
to mix the separated fuel and oxygen. 

Figures 29 and 30 show the previously mentioned variables at 20° ATDC. 

The piston drives the flow throughout most of the cylinder, and the injection 
vortex is somewhat diminished. Near the cylinder wall a new region of combus- 
tion occurs (near the middle of the plot). Figure 30 shows a localized 
high-temperature region. 

Figures 31 and 32 show the flow patterns at 30° ATDC. Figure 31 shows a 
nearly piston-driven flow with the remains of the injection vortex. The cylin- 
der head side of the vortex carries small amounts of fuel into the oxygen-rich 
outer regions. A localized temperature spike in figure 32 shows the presence 
of combustion. Gradually fuel diffuses into the middle regions of the cylin- 
der that were once occupied solely by combustion products. 

Figures 33 and 34 show flow patterns at 40° ATDC, similar to those 
described at 30° ATDC. The remains of the injection vortex carry fuel along 
the cylinder head into the oxygen-rich region causing localized combustion 
regions. 

Figures 35 to 40 show the previously mentioned variables at crankshaft 
angles of 60°, 80°, and 108°. Figures 35, 37, and 39 show a gradual increase 
in the degree the piston dominates the flow. The remains of the injection vor- 
tex still move fuel along the cylinder head into the oxygen-rich outer re- 
gions. The temperature plots (figs. 36, 38, and 40) show the presence of 
small-scale combustion along the cylinder head. Gradually the flame front 
moves radially with the injection vortex. At 108° ATDC, the flow is nearly 
unidirectional. No large-scale motions present in the cylinder are 
capable of completely mixing the remaining fuel and oxygen. 

Approximately the last 88° (20° to 108° ATDC) show very similar behav- 
ior. Gradually the piston begins to drive the flow in a unidirectional man- 
ner. Only one vortex persists, the injection vortex, and it gradually moves 
fuel from the fuel-rich center region toward the oxygen-rich outer cylinder 
along the cylinder wall. The temperature in the cylinder slowly decreases due 
to adiabatic expansion except for localized regions of combustion near the cyl- 
inder head. Gradually thermal boundary layers form as heat is lost to the cy 1 — 
i nder wal 1 s . 

Figures 41 to 45 show plots of the turbulent kinetic energy at 350°, 

360°, 10°, 40°, and 108° ATDC, respectively. Figure 41 shows the turbulence 
levels before the beginning of combustion. Nearly uniform levels of turbu- 
lence exist throughout the cylinder except near the injector, where the injec- 
tion process generates levels of turbulence an order of magnitude larger than 
the background levels. 

Figure 42 shows the turbulence levels immediately after the explosion 
premixed combustion period at 360° ATDC. As in figure 42, the presence of 
three vortices is clearly seen as the local levels of turbulence increase in 
the shear layers associated with the three vortices. These levels are an 
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order of magnitude larger than the background levels. Figure 43 shows the tur- 
bulence at 10° ATDC. The injection vortex and the oxygen-rich-region vortex 
still are strong enough to generate turbulence but global levels decay. 

Figures 44 and 45 show the turbulence levels at 40° and 108° ATDC, respec- 
tively. No substantial large-scale motion generates turbulence in the cylin- 
der, and levels are nearly uniform throughout the cylinder. Only near the 
walls can turbulence generation be observed. 

Figure 46 shows the dissipation rate of turbulent kinetic energy at 350° 
ATDC. Inspection of this figure indicates that the dissipation rate mimics 
the turbulent kinetic energy level. Large levels of turbulence produce larqe 
dissipation rates. 


CONCLUSIONS 

Calculations of the firing two-stroke diesel engine for a 22.5° swirl 
angle were reported, and some conclusions can be summarized as follows: 

1. During the initial period of the combustion cycle (ignition delay 
period), hydrodynamics controls the flow which is primarily driven by the 

vor ^ ex- Large regions of partially premixed fuel are formed around 
the liquid spray region. Vaporization cooling keeps combustion from immedi- 
ately occurring in these regions. 

2. Explosive premixed combustion occurs throughout a large portion of the 
cylinder once combustible temperatures are achieved in the cylinder due to com- 
pression of the gases by the piston and the local fuel/air ratio are within 
the flammability limits. This period lasts from 8° to 10° ATDC. 

3. Large-scale pressure gradients drive the flow immediately after the 
explosive premixed combustion period. The large flow velocities quickly dis- 
perse the remaining droplets throughout the cylinder. 

4. Large temperature rises in the cylinder quickly vaporize a large per- 
centage of the droplet mass which creates (in this geometry) large fuel-rich 
regions near the cylinder centerline. These fuel-rich regions slowly move and 
diffuse toward the oxygen-rich region located near the cylinder wall. 

5. After the explosive period, a diffusion-controlled combustion period 
begins in which Injected fuel drives the flow near the Injector and fuel and 
oxygen move and diffuse together as a result of large-scale fluid motion. 
Diffusion-controlled combustion continues to occur in localized regions 
throughout the cylinder for the remainder of the expansion stroke. 
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TABLE. I. - ENGINE VARIABLES INVESTIGATED 



AXIAL, x 

FIGURE 1. - GEOMETRY OF DIESEL ENGINE CHAMBER. (ALL DIMENSIONS ARE IN INCHES.) 


7 














0.00507 
7 156.270 


\ tt < * << ^^:^;| 




******** * * ♦ < 


*♦♦♦** 


mi 

* * i 

Hit* ^ ^*0*1 

^ ^ ^ ^ ^ ^» ^ <F- <y & tt $ 


^ ^ ) 4. « » <}. , 

71 I* * * U ^ 'l ^ O^ <)- , 


m* ,* u 7 ^ ^ ♦ <► 4 a * * 

"* M * M * ^ ^ * ♦ ♦ 7- 4 A 7 * 

77V?^? ^ ^ ^ ^ ^ ^ /) il ^ ^ f 

*7 7 77 A A £ £ £ £ ^ 

iinii^rr u i t mtl 

* ** * * « 7 r r fM 

777777 7 7 7777 * 

* 7 7 7 A 7 7 7777* 

74747 7 7 4777 # 

77777 7 7 77777 



uj cn 
*- \ 


uj i— 

_j — • 
CL Z 
O CD 
CX C 
Q X 


0.000 
0.00000 

FIGURE 2. - GAS PHASE VELOCITY VECTOR PROFILES AND DROPLET 
LOCATIONS AT 350° ATDC FOR 22.5° SWIRL ANGLE. 
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FIGURE 6. - GAS PHASE VELOCITY VECTOR PROFILE AND DROPLET 
LOCATIONS AT 355° ATDC FOR 22.5° SWIRL ANGLE. 
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FIGURE 8. - GAS PHASE VELOCITY VECTOR PROFILES AND DROPLET 
LOCATIONS AT 356° ATDC FOR 22.5° SWIRL ANGLE. 
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FIGURE 10. - GAS PHASE VELOCITY VECTOR PROFILES AND DROPLET 
LOCATIONS AT 357° ATDC FOR 22.5° SWIRL ANGLE. 
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FIGURE 12. - GAS PHASE VELOCITY VECTOR PROFILES AND DROPLET 
LOCATIONS AT 358° AT DC FOR 22.5° SWIRL ANGLE. 
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FIGURE 13. - TEMPERATURE IS0C0NT0URS AT 358° AT DC FOR 22.5° 
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FIGURE 14. - GAS PHASE VELOCITY VECTOR PROFILES AND DROPLET 
LOCATIONS AT 359° ATDC FOR 22.5° SWIRL ANGLE. 
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FIGURE 16. - 6AS PHASE VELOCITY VECTOR PROFILES AND DROPLET 
LOCATIONS AT 360° ATDC FOR 22.5° SWIRL AN6LE. 
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FIGURE 17. - TEMPERATURE ISOCONTOURS AT 360° ATDC FOR 22.5° 
SWIRL ANGLE. 
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360° ATDC FOR 22.5° SWIRL ANGLE. 
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FIGURE 22 . - GAS PHASE VELOCITY VECTOR PROFILES AND DROP- 
LET LOCATIONS AT 4 ° ATDC FOR 22 . 5 ° SWIRL ANGLE. 
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SWIRL ANGLE. 
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LET LOCATIONS AT 15° ATDC FOR 22.5° SWIRL ANGLE. 
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FIGURE 28. - TEMPERATURE 1S0C0NT0URS AT 15° ATDC FOR 22.5° 
SWIRL ANGLE. 
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FIGURE 29. - GAS PHASE VELOCITY VECTOR PROLILES AND DROPLET 
LOCATIONS AT 20° ATDC FOR 22.5° SWIRL ANGLE. 
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FIGURE 30. - TEMPERATURE ISOCONTOURS AT 20 u ATDC FOR 22. 5 U 
SWIRL ANGLE. 
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FIGURE 32. - TEMPERATURE ISOCONTOURS AT 30° AT DC FOR 22.5° 
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FIGURE 33. - GAS PHASE VELOCITY VECTOR PROFILES AT 40° ATDC 
FOR 22.5° SWIRL ANGLE. 
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FIGURE 34. - TEMPERATURE IS0C0NT0URS AT 40° ATDC FOR 22.5° 
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FIGURE 35. - GAS PHASE VELOCITY VECTOR PROFILES AT GO 0 ATDC 
FOR 22.5° SWIRL ANGLE. 
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FIGURE 36. - TEMPERATURE ISOCONTOURS AT 60° AT DC FOR 22.5° 
SWIRL ANGLE. 
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FOR 22.5° 


p p p p p 

p p p p p 

p p p p p 

p p p p p 

p p p p p 

p p p p p 

p p p p p 

p p p p p 

p p p p p 


- \ 


a lu 
o 

t- 3 
LU l— 
_J — 
Q. Z 
O CD 

a c 

a 3: 


0.000 
0.00000 

- GAS PHASE VELOCITY VECTOR PROFILES AT 80° ATDC 
SWIRL ANGLE. 


0.00507 

7156.270 



8 7 8 2 . 8 7 9 


4-4-6 7°. 4-9 6 


2272.4-3 1 


1 1 5 5*. 8 9 3 

LU 

O' 

3 

£ 5 8 7*. 9 5 5 

LU 

CL 

£ 

UJ 

*" ■ 

419.332 


FIGURE 38. - TEMPERATURE IS0C0NT0URS AT 80° ATDC FOR 22.5° 
SWIRL ANGLE. 
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FIGURE 39. - GAS PHASE VELOCITY VECTOR PROFILES 
FOR 22.5° SWIRL ANGLE. 
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FIGURE 41. - NATURAL LOGARITHM OF TURBULENT KINETIC ENERGY 
PROFILES AT 350° ATDC FOR 22.5° SWIRL ANGLE. 
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FIGURE 42. - NATURAL LOGARITHM OF TURBULENT KINETIC ENERGY 
ENERGY PROFILES AT 360° ATDC FOR 22.5° SWIRL ANGLE. 
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FIGURE 43. - NATURAL LOGARITHM OF TURBULENT KINETIC ENERGY 
PROFILES AT 10° ATDC FOR 22.5° SWIRL ANGLE. 
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FIGURE 44. - NATURAL LOGARITHM Of TURBULENT KINETIC ENERGY 
PROFILES AT 40° ATDC FOR 22.5° SWIRL ANGLE. 


FIGURE 45. - NATURAL LOGARITHM OF TURBULENT KINETIC ENERGY 
PROFILES AT 108° ATDC FOR 22.5° SWIRL ANGLE. 
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FIGURE 46. - NATURAL LOGARITHM OF DISSIPATION RATE OF TUR- 
BULENCE KINETIC ENERGY AT 350° ATDC FOR 22.5° SWIRL ANGLE. 
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